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ABSTRACT 

I present results from an optical spectroscopic investigation of the binary 
system HD 42401 (V1388 Ori; B2.5 IV-V + B3 V). A combined analysis of 
y— band photometry and radial velocities indicates that the system has an orbital 
period of 2.18706 ± 0.00005 days and an inclination of 75.5 ± 0.2 degrees. This 
solution yields masses and radii of Mi = 7.42 ± 0.08 M© and Ri = 5.60 ± 0.04 
Rq for the primary and M2 = 5.16 ± 0.03 Mq and R2 = 3.76 ± 0.03 Rq for the 
secondary. Based on the position of the two stars plotted on a theoretical H-R 
diagram, I find that the age of the system is > 25 Myr and that both stars appear 
overluminous for their masses compared to single star evolutionary tracks. A fit 
of the spectral energy distribution based on photometry from the literature yields 
a distance to HD 42401 of 832 ± 89 parsecs. 

Subject headings: binaries: spectroscopic - stars: early-type - stars: fundamental 
parameters - stars: individual (HD 42401) 



1. Introduction 

Two fundamental parameters in stellar astrophysics are the masses and radii of stars. 
Especially lacking are accurate (< 2%) determinations of these quantities for O- and B- 
type systems. Data of this kind are important for tests of stellar models, with far-reaching 
implications such as the modeling of stellar components of galaxies. Here I describe mea- 
surements of these parameters for the B-star binary HD 42401. This effort is part of an 
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ongoing program to determine the masses and radii of O- and B -type stars. Previou sly, I 
obtained parameters for the LMC 0-star binary [L72] LH 54-425 (jWiUiams et aLll2008l ) and 
will continue with more analyses of eclipsing spectroscopic binary systems containing O- and 
B-type stars. 



The star HD 42401 was first classified as B2 V by IWalborru (jl97ll ). Since then, it 
has b een used as a spectral standard for the B2 V type stars in a numbe r of publica- 
tions JWalborn fc Fitzpatricklll99ol : Ide Mello et al.l[200o[ iBagnuolo et allboOlL The eclips- 
ing nature of its light curve was first noticed with the Hipparcos satellite (HIP 29321), 
where it was classified as an Algol-type eclipsing binary and given the mor iiker V1388 Ori 
( Kazarovets et al. I ll999h . Three radial velocity measurements of HD 42401 by lFehrenbach et al, 
( 119971 ) show a range from -103 to 36 km s~^, but no further investigations were made into 
this probable velocity variable. 

Presented here are the analyses of radial velocities from optical spectra obtained in 
2008 January (§2). Section 3 contains a description of the measurement of radial velocities 
from the data. In §4 I describe tomographic reconstructions of the individual component 
spectra. Covered in §5 is the combined light curve and radial velocity curve solution. In 
§6 I finish with a discussion of the distance to the system, its fundamental parameters, and 
other consequences of the analysis. 



2. Observations 



I obtained 29 spectra of HD 42401 with the Kitt Peak National Observatory (KPNO) 
0.9-m coude feed telescope during an observing run from 2008 Jan 09 to 2008 Jan 18. These 
spectra were obtained with the long collimator and grating A (632 grooves mm~^ with a blaze 
wavelength of 6000 A in second order with order sorting filter 4-96) . The first two nights of 
observations were made with the TIKB detector (1024x1024 pixel array with 24 square /im 
pixels) and resulted in a resolving power of R = X/6X = 13,500 with wavelength coverage 
from 4315-4490 A. The rest of the observing run made use of the F3KB detector (3072x1024 
pixel array with 15 square fim pixels) that resulted in a lower resolving power of R = X/6X = 
11,500 with wavelength coverage from 4250-4570 A. The wavelength coverage was chosen 
to include the He I A4471 and Mg ll A4481 lines that are good temperature indicators for 
B-type stars plus the H7 line that is sensitive to gravity (linear stark effect). Exposure times 
were 1200 s giving a S/N ~ 100 pixel"^. Numerous comparison spectra were obtained for 
wavelength calibration, and many bias and flat field spectra were also obtained each night. 
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The spectra were extracted and calibrated using standard routines in and then each 

continuum rectified spectrum was transformed to a common hehocentric wavelength grid in 
log A increments. 



Radial Velocities 



The narrow range of wavelength coverage in our spectra limits the lines that may be 
used for radial velocity analysis. I measured radial velociti es from three liri es, He I A4387, 
He I A4471, and Mg ll A4481 via a template-fitting scheme (IGies et al.l 120021 ) that measures 
velocities by using model templates weighted by a fiux ratio to match both the shifts and line 
depths in the observed spectra. There is no evidence of emission or intrinsic line asymmetries 
in these lines or H7 A4340. 



The BSTAR2006 grid of stellar models from lLanz fc Hubenyl (120071 ) was used to derive 
template spectra. These models are based on the line blanketed, non-LTE, plane- parallel, 
hydrostatic atmosphere cod e TLUSTY and the r adiative transfer code SYNSPEC ( iHubeny 
19881 : iHubeny &: Lanzlll995l : iHubeny et al.lll998l ). In finding templates, initial values were 
used for temperatures, gravities, projected rotational velocities, and fiux contributions from 
each star. These parameters for model templates were then checked by comparing the three 
lines used in radial velocity measurements against the tomographically reconstructed spectra 
of each star (^. The parameters were changed and new templates made after initial fits 
to the light curve and radial velocity curves (§5) indicated slightly different values were 
more appropriate. The velocity analysis was then performed again until the best fit was 
obtained. In this way, I was able to obtain an estimate of the monochromatic fiux ratio in 
the blue spectra of F2/F1 = 0.25 ± 0.05 based upon the relative line depths of the spectral 
components. 

The template fitting scheme also needs preliminary estimates for the velocities of each 
component. To obtain these, spectra which clearly showed two sets of spectral lines were an- 
alyzed. Midpoints of spectral features were used to make crude estimates of radial velocities 
for each star. Preliminary orbital parameters were obta ined by using the nonlinear, least- 
squares fitting program of iMorbey fc Brosterhud (119741 ). The preliminary velocities from 
this initial orbital solution were starting points for performing a nonlinear, least-squares fit 
of the composite profiles with the template spectra and calculating the shifts for each star. 
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The fitting scheme searches a region of ±10 A around the rest wavelength of each line, so 
preliminary velocities are mere starting points and need not be highly accurate. The three 
values for radial velocity from each spectrum were averaged, and the standard deviation of 
the mean value was calculated. Each of these values are listed in Table [T] for the primary 
and secondary stars. Also listed are the orbital phase for each observation and the observed 
minus calculated values for each data point. Phase zero is defined as the time of inferior 
conjunction of the primary star, Tic,i (time of secondary minimum in the light curve). Figure 
[T] shows two sample spectra at times of near quadrature according to this ephemeris. 

During eclipse phases, velocity measurements may be affected by the Rossiter-McLaughlin 
effect in which the center of light of the eclipsed star will appear redshifted on ingress and 
blue shifted on egress because of the rotational Doppler shifts of the visible portions. This ef- 
fect will result in narrower line profiles as the red and blue shifted regions of the eclipsed star 
are blocked, and it is more prominent in broad- lined, rapidly rotating systems. HD 42401 
is a short period binary having relatively narrow lines, therefore the effect should be small. 
However, since the template method does not account for eclipses, deviations from orbital 
motion near eclipse phases are expected. 



4. Tomographic Reconstruction 



The Doppler tomography algorithm of iBagnuolo et al.l (119941 ) was used to separate the 



primary and secondary spectra of HD 42401 for the F3KB spectra. This iterative method 
uses the 24 observed composite spectra, their velocity shifts, and an assumed monochromatic 
flux ratio to derive individual component spectra. The best flux ratio was the one that best 
matched line depths in the reconstructions with those in the model spectra. Figure [2] shows 
the reconstructed spectra for the primary and secondary with the best fit models over- 
plotted. The relative depths of He I A4471 and Mg ll A4481 are good temperature indicators 
throughout the B-star sequence. Specifically, the He I A4471 line gets weaker while the 
Mg II A4481 line gets stronger as temperature decreases, as is seen in the spectrum of the 
secondary compared to that of the primary. The set of spectra from the F3KB instrument 
has only one spectrum that is close to an eclipse. Because there are many more spectra 
outside eclipse used in the tomographic reconstruction, the final reconstructed spectra are 
insensitive to details of the Rossiter-McLaughlin effect. Indeed, tomographic reconstructions 
with the one eclipse spectrum omitted are negligibly different from those presented in Fig. 

m 



The final reconstructed spectra were fit with TLUSTY/SYNSPEC model synthesis spec- 
tra (see IJ3]). Fits to all three lines made in the velocity analysis were used to estimate the 
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projected rotational velocity Vsini, and the temperature and surface gravity of each star 
were estimated by c omparing the reconst ructed and model profiles for a grid of test values 
(details are given in lWilliams et al.ll2008l ). These parameter values are listed in Table [21 



The narrow wavelength range of the spectra is insufficient to attempt classical spectral 
typing. However, I can arrive at an estimate of the spectral types for each star in the 
HD 42401 system by comparing the derived effective temperat ures with a sp e ctral type 
versus effective temperature relation. According to Table 2 of iBohm-Vitensd (jl98ll ) the 
effective temperature and gravity (§5) of the primary of HD 42401 is most consistent with 
a B2.5 IV-V star while the secondary matches most closely with a B3 V star, and these 
classifications are listed in Table [2j 



5. Combined Radial Velocity and Light Curve Solution 



The All Sky Automated Survey (ASAS-3; iPojmanskil |2002| ) database contains light 
curves for ~39,000 previously unknown variable stars. I extracted the l^-band light curve for 
HD 42401 (ASAS 061059+1159.7) from this catalog, removing points deemed lower quality 
by the data reduction pipeline used by the ASAS. 



I used the Eclipsing Light Curve (ELC) code (lOrosz &: Hauschildtll2000l ) to find orbital 
and astrophysical parameters for the HD 42401 system. ELC fits the radial velocity and light 
curves simultaneously, giving a joint orbital ephemeris based on both sets of data. The best 
fit for the radial velocity curve is shown in Figure [3]^ upper panel) along with the observed 
minus calculated residuals (lower panel). ELC treats calculation of v e lociti es during eclipse 
via a fiux- weighted velocity centroid as described in IWilson fc Sofia! (119761 ). The Rossiter- 
McLaughlin effect is clearly seen in the best fit as a slight redshift going into eclipse and a 
slight blueshift leaving eclipse and the model velocity curves match the data well. The best 
fit for the l^-band light curve is shown in Figure HI Figures [5] and [6] show the details of the 
fit for the light curve in the region of the primary and secondary eclipses, respectively. 

The genetic optimizer mode of ELC was used initially to explore wide ranges of values 
for the period, epoch of inferior conjunction of the primary Tic,i, inclination, mass ratio, 
primary velocity semiampl itude, and Roche lobe fillin g factor for each star. The Roche lobe 
filling factor is defined by lOrosz fc Hauschildtl (120001 ) as the ratio of the radius of the star 
toward the inner Lagrangian point (Li) to the distance to Li from the center of the star, 

f — '^point/'^Ll • 

The use of ELC included fixing the temperature of each star to the values found in 
the tomographic reconstructions (§!]) and also fixing the radius ratio {R1/R2) based upon 
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the temperatures, surface fluxes, and monochromatic flux ratio of the two stars. Non-zero 
eccentricities for the HD 42401 system were explored during fitting but rejected based on 
the higher values for those fits. 

To estimate the uncertainties based on our best fit, the values of the seven fitted pa- 
rameters were varied in the calculation of ~ 3 x 10^ light and radial velocity curves. The 
well-explored surface was then projected as a function of each fitted parameter or astro- 
physical parameter of interest. The lowest value is found for each parameter, and the 
l-cr uncertainty may be estimated by the region where ^ Xmin + 1- These values and 
uncertainties are listed in Table [3] for the orbital parameters and Table H] for the astrophysi- 
cal parameters. Also listed are i?efr, the effective radius of a sphere with the same volume, 
i?poie, the polar radius of each star, and -Rpoint, the radius of the each star toward the inner 
Lagrangian point. Our derived period of P= 2.18706 ± 0.00005 days for HD 42401 agree s 



well with the value found by the Hipparcos satellite of 2.18709 ± 0.00050 days flESAlll997h . 



6. Discussion and Conclusions 

The goal with this work was to determine the masses and radii of the stars in the 
HD 42401 system to great accuracy. I have determined the masses of the two stars to 
1.1% and 0.6% and the radii of the stars to 0.7% and 0.8%, for the primary and secondary 
respectively. Armed with these numbers and the results from Table HJ I can consider the 
details of the HD 42401 system. 

As is seen in Table HJ both stars are well within their Roche radii but experience tidal 
distortion that is evident in the light curve (Fig. Sj). The rotational velocities derived from the 
tomographic reconstructions (Table [2]) match very well with the synchronous rotation values 
found by ELC (Table H]), indicating that this system has achieved synchronous rotation, and 
is therefore not very young. 

To make an estimate of the age of the system, I used the effective temperatures from 
Table[2]and radii from Table|l]to plot the two stars of HD 42401 on a theoretical H-R diagram 
and compare their locations to evolutionary tracks. The res ult is shown in Figure [TJ plotted 



against evolutio nary tracks for stars of 5, 7, and 9 Mq from ISchaller et al.l (119921 ) as well as 



isochrones from lLejeune fc Schaereii (120011 ) for solar metallicity with ages of 21.9, 25.1, 27.5, 
and 31.6 Myr. The location of the stars is most consistent with an age of ~25 Myr. The 
model tracks shown in Figure [7] are for non-rotating stellar models. The present ratios of 
spin angular velocity to critical angular velocity are Q/Qcrit = 0.45 and 0.31 for the primary 
and secondary, respectively (assuming synchronous rotation), and the evolutionary tracks 
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for such moderat e rotation rates are o nly slightly steeper and more extended in time than 
those illustrated (lEkstrom et al.l l2008l ) . Thus, our derived age may slightly underestimate 
the actual value. 

Both stars appear (in Fig. [7]) to be overluminous f or the derived ma sses of 7.42 Mq for 
the primary and 5.16 Mq for the secondary. Table 3 of iHarmaned (119881 ) lists astrophysical 
parameters for stars as a function of spectral type for main sequence stars based on empirical 
data from eclipsing binaries. The mass and effective temperature of the primary fit between 
the listed values for spectral types B2 (mean of 8.6 Mq) and B3 (mean of 6.1 M©), but the 
radius is much larger than the means for comparable spectral types and matches a BO. 5 star 
(mean of 5.5 Rq). The mass of the secondary is consistent with the B4 spectral type (mean 
of 5.1 Mq) while the effective temperature and radius appear more consiste nt with the B3 



spectr al type. In a study of eclipsing binaries in the Small Magellanic Cloud, iHilditch et al. 



(120051 ) found several systems of comparable mass that, like HD 42401, are overluminous 
compared to model predictio ns. However, our findings for HD 42401 seem to conflict with 
the results of iMalkovl (120031 ) who shows that early B-type stars that are in close systems 
and rotate more slowly t han sing l e star s are on average smaller than those same single stars. 
In a subsequent paper, iMalkovl (120071 ) studied well-separated binaries in an effort to use 
the properties of their component stars for a more direct comparison with single stars. His 
resulting mass- luminosity-radius relations, when applied to our results for HD 42401, predict 
less luminous, hotter and smaller components. This is perhaps not surprising, due to the 
age of HD 42401 and the evolution of its components from the zero age main sequence. 

I can also estimate the distance to the system by fitting a spectral energy distribution 
(SED) to various photometric measurements. The spectra were not flux calibrated, so to 
create an SED, I must rely on photometry performed on the system. HD 42401 is a relatively 
bright {y ~ 7.4) system and has thus been well studied. The problem arises in estimating 
the phase at which a particular observation in the historical literature was made. Two obser- 
vations of HD 42401 were taken by lUE at = 0.12 ± 0.07 based on our orbital ephemeris. 
The uncertainty in this estimate comes from the uncertainties in our values for Tic,i and 
period and the number of orbits between the lUE observations and our observations. The 
range in uncertainty for the orbital phase of the lUE data brings it close to eclipse, and this 
was supported by initial SED fits including the lUE data that indicated a lower UV flux 
than indicated by other measurements. There are various other measurements th at also were 



suffic iently close to eclipses to be omitted in SED fitting. Fortunately, 2MASS (ICutri et al. 
2003! ) measurements were obtained at = 0.16±0.03 and are therefore far enough away from 



eclipse, within uncertainties, to be used. The only other points used were Johnson UBV 
magnitudes. To obtain these values, I used the value at quadrature from our light curve of 
V = 7.40 and applied colors of [B — V) = —0.042 and {U — B) = —0.620 which are averages 
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of several observations that are listed by iMermilliod fc MermilliodI (Il994l ). The SED fit is 
show n in Figure [8] along wi th the U,B,V, J, H, Kg magnitudes. I used two model spectra 
from iLanz fc Hubenyl (120071 ) matching our values in Table [2] for the effective temperatures 
and gravities of each star, and scaled in the blue by the flux ratio of 0.25 ± 0.05 found from 
the tomographic reconstructions. This fit of the SED results in a limb darkened, angular 
diameter for the primary of ^ld = 62.7±1.1 /las with a reddening of E{B — V) = 0.15±0.01 
mag and a ratio of total-to-selective extinction of R = 3.24 ± 0.10. By directly comparing 
this angular diameter with the value for the radius of the primary, I estimate the distance of 
the system to be d = 832 ± 89 pc T his distance and reddening are in excellent agreement 
with the values in ISowen et aD J2008h of E{B -V) = 0.15± 0.02 mag and d = 0.8 kpc. 



HD 42401 has Galactic coordinates of £ = 197?64 and b = -3?33 flReedl 120051 ). This 
is close t o th e Gala ctic open cluster NGC 2169 at i = 195?61 and b = — 2?93 (separation 
~ 2?07). lAbtl (119771) fourid that the earliest spectral type of the cluster members of NGC 2169 
was B2 III. Ijeffries et all koO'^ find a reddening value of E{B -V) = 0.20 ±0.01 mag and a 
distance of ~ 1060 pc to NGC 2169. Proper motions for the objects are similar as well, with 
HD 4 2401 having jj^ cos 6 = -2.25 ± 1.41 mas yr"^ and fis = -2.71 ± 0.53 mas yr"^ JeSA 



1997D and NGC 2169 having /x„ cos 6 = -2.17 ± 1.41 mas yr^^ and fis = -2.52 ± 0.53 mas 
yr~^ (IKharchenko et al.ll2005l ). The systemic velocity deriv ed here for HD 42401 is ^ 14 km 
s~^ and the average of 8 stars in NGC 2169 is 22 km s~^ (IKharchenko et al.ll2005l ). While 
the proximity of HD 42401 to NGC 2169 is intriguing and the fact that the earliest star has a 
similar spectral type to the components of HD 42401, actual cluster membership is unlikely. 
The age of NGC 2169 in the literature ranges from 9 ±2 Myr Jjeffries et al1l2007[ ) to 7.8 Myr 
(IKharchenko et al.l 120051). too voung to comfort ably fit with the data for the components of 
HD 42401 in Figure [71 IKharchenko et al.l (120051 ) determined the corona radius of NGC 2169 
by matching stellar densities of the cluster to that of the background and found an angular 
size of 0.15 degrees. The distance on the sky of HD 42401 from NGC 2169 is about 14 times 
this radius. The difference in distances to the two objects and their separation on the sky 
make cluster membership even more unlikely. However, it is possible that a wave of star 
formation associated with the Orion arm of the Galaxy swept through the region, leaving 
HD 42401 behind as it moved on and eventually led to the creation of NGC 2169. 
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Table 1. HD 42401 Radial Velocity Measurements 



Date 
(HJD-2,400,000) 


Orbital 
Phase 


Vi 

(km s-i) 


(km s-l) 


(O - C)i 
(km s-i) 


(km s-i) 


(72 

(km s-l) 


(O - C)2 

(km s-i) 


54474.738 


0.548 


-36.3 


5.4 


6.9 


90.6 


3.8 


12.9 


54474.836 


0.593 


-71.7 


3.3 


-1.3 


138.1 


13.8 


4.9 


54476.690 


0.441 


82.7 


0.4 


2.3 


-45.7 


6.3 


20.4 


54476.735 


0.461 


65.6 


3.3 


-1.6 


-21.1 


8.3 


18.1 


54476.775 


0.480 


53.8 


1.6 


5.7 


-35.8 


21.7 


-20.9 


54477.672 


0.890 


-84.6 


1.8 


-3.8 


152.6 


6.9 


1.2 


54477.720 


0.912 


-70.4 


4.0 


-6.5 


129.4 


4.7 


2.4 


54477.800 


0.949 


-41.4 


3.0 


-8.9 


84.4 


2.0 


-13.8 


54478.660 


0.342 


139.3 


0.9 


-3.9 


-160.9 


5.4 


8.2 


54478.676 


0.349 


133.3 


2.0 


-6.3 


-160.1 


5.6 


3.6 


54478.718 


0.368 


125.0 


1.9 


-3.3 


-156.2 


5.4 


-9.2 


54478.779 


0.396 


109.6 


4.1 


0.3 


-114.4 


11.0 


4.9 


54478.840 


0.424 


87.6 


2.0 


-1.6 


-76.1 


16.9 


10.9 


54479.667 


0.802 


-128.4 


0.9 


-1.9 


218.7 


4.6 


0.6 


54479.718 


0.825 


-115.9 


1.2 


2.3 


214.2 


5.3 


8.3 


54479.760 


0.845 


-109.8 


0.9 


-0.8 


194.8 


2.0 


2.4 


54479.825 


0.874 


-90.2 


0.5 


1.3 


174.5 


2.6 


7.5 


54480.635 


0.245 


163.8 


0.7 


-1.5 


-197.1 


6.8 


6.6 


54480.684 


0.267 


163.2 


3.8 


-1.6 


-196.5 


5.7 


6.2 


54480.755 


0.300 


155.8 


2.3 


-3.3 


-194.5 


4.2 


-1.0 


54480.802 


0.321 


151.8 


2.0 


-0.5 


-186.8 


4.3 


-3.9 


54481.635 


0.702 


-127.3 


2.5 


1.6 


240.8 


6.7 


20.8 


54481.680 


0.722 


-135.7 


2.4 


-2.7 


226.5 


6.5 


0.0 


54481.728 


0.744 


-135.8 


0.6 


-1.1 


222.2 


2.5 


-7.4 


54481.794 


0.775 


-135.3 


2.7 


-2.5 


234.0 


4.6 


6.9 


54483.646 


0.622 


-94.4 


2.2 


-3.3 


171.4 


5.9 


8.2 


54483.687 


0.640 


-108.5 


1.9 


-5.8 


180.2 


2.5 


-0.1 


54483.757 


0.672 


-120.3 


2.0 


-1.5 


200.8 


2.9 


-3.6 


54483.839 


0.710 


-129.4 


2.0 


1.3 


223.4 


3.6 


0.5 



Table 2. Tomographic Spectral Reconstruction Parameters 



Parameter Primary Secondary 



Spectral Type . . . 


B2.5 IV- 


B3 V^ 


Teff (kK) 


20.5 ± 0.5 


18.5 ± 0.5 


log g (cgs) 


3.75 ± 0.25 


4.00 ± 0.25 


V sin i (km s~^) 


125 ± 10 


75 ± 15 


F2/F1 (blue) 


0.25 ± 0.05 



^These spectral types are estimated from de- 
rived values of Tes and log g. 
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Table 3. Circular Orbital Solution for HD 42401 



Element 



Value 



P (days) 2.18706 ± 0.00005 

Tic,i (HJD-2,400,000) . . . 54477.9130 ± 0.0002 

Ki (km s-^) 151.4 ± 0.3 

X2(kms-i) 217.9 ±1.0 

71 (km s-^) 15.32 ± 0.07 

72 (km s"^) 12.9 ± 0.2 

rms (primary) (km s~^) . 3.7 

rms (secondary) (km s~^) 9.7 

rms (photometry) (mag) . 0.007 



Table 4. ELC Model Parameters for HD 42401 



Parameter 



Primary 



Secondary 



Inclination (deg) . . . 




75.5 ± 0.2 




M (Mo) 


7.42 


± 0.08 


5.16 


± 0.03 


-Reflf (-R©) 


5.60 


± 0.04 


3.76 


± 0.03 


-Rpole'^ {Rq) 


5.40 


± 0.02 


3.70 


± 0.01 


-Rpoint*^ (-^0) 


5.97 


± 0.02 


3.87 


± 0.01 


Vsync sin i (km s~^) 


125.5 


± 0.9 


84.3 


± 0.6 


log g (cgs) 


3.812 


± 0.005 


3.999 


± 0.006 


Filling factor 


0.674 


± 0.005 


0.588 


± 0.010 






16.48 


± 0.05 




F2/F1 (blue) 




0.27 ± 0.04 





^Polar radius. 

^Radius toward the inner Lagrangian point. 
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Fig. 1. — Two spectra of HD 42401 obtained near quadrature phases (offset for clarity). 
Solid lines indicate the positions of lines of the primary component and dashed lines indicate 
those for the secondary component. The lines marked are those used in the velocity analysis. 
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Fig. 2. — Tomographic reconstructions of the components of HD 42401 based on 24 spectra 
obtained with the F3KB instrument during January 2008 at KPNO. The top sohd hne rep- 
resents the reconstructed spectrum of the primary, the bottom sohd hne is the reconstructed 
secondary spectrum offset by 0.3 for clarity. Over-plotted for both are the model spectra for 
each (dashed lines). The stellar parameters for the model spectra are given in Table El 
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Orbital Phase 



Fig. 3. — Radial velocity curves for HD 42401. Primary radial velocities are shown by filled 
dots and secondary radial velocities are shown by filled triangles with associated uncertainties 
shown as line segments for both. The solid line is the best fit solution for the primary and 
the dashed line is the same for the secondary. The lower panel shows the observed minus 
calculated values for each measurement with uncertainties. 
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Fig. 4. — l^-band light cur ve for HP 42401 . These data were taken from the All Sky 
Automated Survey database (1Pojmanskill2002l ) and are presented here in phase according to 
our best fit solution. The model is the thick gray line and the data are the filled dots with V 
uncertainties represented by line segments. Phase zero corresponds to inferior conjunction 
of the primary star. 
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Fig. 5. — y-band light curve for HD 42401 around the time of primary echpse. As with Fig. 
m the model is the thick gray line and the data are filled dots with line segments representing 
uncertainties. Residuals are distributed as expected around the best fit. 
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Fig. 6. — V^-band light curve for HD 42401 around the time of secondary echpse. As with Fig. 
m the model is the thick gray line and the data are filled dots with line segments representing 
uncertainties. Here, as in Fig. [5], the residuals are distributed as expected around the best 
fit. 
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Fig. 7. — A theoretical H-R diagram showing the location of the primary star [filled circle) 
and secondary star [filled triangle) of HD 42401 including uncerta inty regions for each. Also 
plotted are evolutionary tracks for sta rs of various masses from ISchaller et al.l (Il992l ) and 
isochrones [vertical dashed lines) from lLejeune &: Schaererl (120011 ) for solar metallicity with 
ages of 21.9, 25.1, 27.5, and 31.6 Myr going from left to right. The positions of the two stars 
are consistent with an age of ~25 Myr. 




Fig. 8. — The spectral energy distribution and fit for the combined hght of the HD 42401 
components {solid line) to Johnson U,B,V, J, H, Kg photometry {open circles). 



